Irreversibility in response to forces acting on graphene sheets 
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The amount of rippling in graphene sheets is related to the interactions with the substrate or 
with the suspending structure. Here, we report on an irreversibility in the response to forces that 
act on suspended graphene sheets. This may explain why one always observes a ripple structure on 
suspended graphene. We show that a compression-relaxation mechanism produces static ripples on 
graphene sheets and determine a peculiar temperature Tc, such that for T < Tc the free-energy of 
the rippled graphene is smaller than that of roughened graphene. We also show that Tc depends on 
the structural parameters and increases with increasing sample size. 
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Introduction — Two-dimensional graphene crystals ^J 
have attracted considerable attention both experimen- 
tally and theoretically, due to their unusual electronic 
properties [2] . Ripples or undulations were first observed 
in freely-suspended graphene flakes in experiments [3]. 
The ripples affect the electronic properties [¥|, such as 
the conductivity and the quantum transport properties. 
Charge inhomogeneities due to the ripples have also been 
observed in graphene [5]. A new method to fabricate 
periodically rippled graphene on Ru(OOOl) under ultra- 
high vacuum conditions was reported in [5]. The ripples 
in graphene were studied theoretically [7], and it was 
claimed that the ripples may be explained as a conse- 
quence of absorbed molecules sitting on random sites. 

There are two points of view on the physics of the rip- 
ples. Meyer et al. [5] proposed that the reproducible ap- 
pearance of the ripples across samples indicates that it is 
an intrinsic effect. They emphasized that the homogene- 
ity and isotropy of the ripples are not compatible with the 
assumption of an incompressible sheet. They estimated 
a local strain of up to 1% for a single-layer flake. Using 
the generic free-energy for the long wavelength deforma- 
tions [8j, Castro Neto and Kim [9], on the other hand, 
argued that graphene may be considered as an atomic 
thin membrane, with its physics being also similar to a 
soft membrane. In this point of view the ripples are not 
intrinsic and can be the results of the environment, such 
as the substrate in the system. 

Recently, direct observation and controlled creation 
of periodic ripples in suspended graphene sheet was re- 
ported jlOj by using both spontaneously and thermally 
generated strains via varying the substrate and annealing 
conditions. The experimental measurements indicated 
that the ripples were induced by the preexisting longitu- 
dinal strains in graphene. It was observed at tempera- 
tures about 500 K that graphene sheets are flat, however, 
upon cooling down to room temperature, the ripples in- 
variably appear. 

In this Letter we report on an irreversibility in re- 
sponse to forces that act on the suspended graphene 



sheets that may explain why one always observes the rip- 
pled structure on graphene in experiments. We show that 
a compression-relaxation mechanism can produce static 
ripples on graphene sheets. We determine a peculiar tem- 
perature Tf, such that for temperatures less than T^ the 
free-energy of the rippled graphene is smaller than that 
of roughened graphene. We also show that for sample 
with size 400 x 200 atoms, Tc ~ 90 K and, moreover, Tc 
is an increasing function of the sample size. 
Theory and Model — We used molecular dynamics (MD) 
simulations with the empirical interatomic interaction 
potential due to Brenner |llj . i.e., carbon-carbon inter- 
action in hydrocarbons that contains three-body inter- 
action. Many-body effects of electron system, on aver- 
age, are considered in the Bernner potential through the 
bond-order term. We employed the Nose-Hoover thermo- 
stat to control the temperature, when we used a canonical 
(NVT) ensemble in the MD simulations. We note that 
although the Brenner potential is not entirely a quantum 
mechanical potential, it predicts the correct mechanical 
properties of the structures with carbon atoms by using 
the classical MD simulations |12j . 

Here, we show that the compression and then re- 
laxation in one or two directions [x (arm-chair) and y 
(zigzag)] of graphene sheet can produce static ripples, 
which means that there is an irreversibility in response 
to forces acting on segments of the graphene sheets. In- 
deed, we have found that if one compress the surface 
in one direction, say the a;— direction (arm-chair), after 
compressing about 0.13%i, where L is the size of the sim- 
ulation sample, the ripples will appear. After the ripples 
emerge, we move back the boundaries to their imimZ posi- 
tions. We then observe that after doing the compression- 
relaxation processes, the ripples survive, hence imply- 
ing that the compression process is not reversible. It is 
worthwhile mentioning that we also simulated tethered 
membranes |l3j and repeated the compression-relaxation 
procedure. We obtained no any indication of irreversibil- 
ity in response to the forces acting on the membranes. 

In the case of graphene, if the compression amount be- 



comes larger than the critical value (here, 0.13%L, which 
also depends on the temperature and the size of simula- 
tion sample), the graphene sheet bends and, therefore, 
no ripple appears. The typical height variance of the rip- 
pled graphene is about 5 A at T = 50 K. Our simulation 
results show that the wavelength of the static ripples do 
not change with the size of the sample. Moreover, we 
observe that the surface roughness (the variance of the 
height fluctuations) does not change after the relaxation 
and, therefore, the ripples are static. Thus, we might 
state that any primary stress on graphene sheet, for ex- 
ample in its preparation in the experiments, can con- 
struct ripples that will survive during the experimental 
measurements (see, for example, |14j). 

Let us first determine the average wavelength of the 
ripples after relaxing the system. To do so we calculate 
the two-dimensional Fourier components of the height- 
height correlations, G(|q|) =< |/i(q)| >. Figure 1 shows 
Gdql) as a function of |q|/goi in logarithmic scales, for 
both the roughened (no ripples) and relaxed states in 
which we have stable ripples. Here qg = 2t:/L, with L 
being the length of graphene in the x— direction. In the 
inset of Fig. 1 the same plot in linear scale for |q| is shown 
to clarify a peak around |q| ~ lOqo that corresponds to 
about 85 A at 50 K. This is the average wavelength of the 
ripples and is near to the value observed experimentally 
[3] and calculated numerically [THITB]. In addition, one 
can derive the wavelength of the ripples by calculating 
the first minimum of the second moments of the height 
increments fluctuations < \h{xi) — h{x2)\'^ > with respect 
to relative distance, \xi — X2I. 

The scale-dependence of < \h{q)\ > is proportional 
to l/|q|", where a ~ 4 at temperature 50 K. Conse- 
quently, our results predict that the bending rigidity term 
prevails with respect to the surface tension in graphene 
at short distances [TTj. Note that the contribution of 
surface tension is a term like T/cr|q| , whereas the con- 
tribution of the bending rigidity is T/K|q| , where a 
and K are the interfacial tension and bending modulus, 
respectively J17j . The exponent a might generally be 
smaller than 4, due to thermal fluctuations, surface ten- 
sion and anharmonic corrections fS^. We estimate k, the 
bending rigidity or bending modulus, using the relation, 
K-i ~ |q|4 < \h{q)\^ > /NkeT. Plotting k vs |q|/go 
shows that the bending rigidity is almost constant for 
20 < |q|/qo < 100, with k ~ 1 eV^i. 

Thus, for a given temperature, the compression- 
relaxation mechanism produces the ripples, and graphene 
has at least two "states" simultaneously, the "rough" or 
normal sheet (no ripple) and the "rippled" structure. The 
question now is, which state is more stable? To answer 
this question one should calculate the free energies of the 
two states and determine which state has a smaller free- 
energy. In what follows we calculate the free-energy dif- 
ference of the rippled and roughened states of graphene 




FIG. 1: (Color online) Log-log plot of < |^t(q)|^ > as a func- 
tion of q/go (log scaled) both in the rough and ripple cases. In 
the inset, < |/!,(q)p > is shown as a function of jqj/go to clar- 
ify a peak around |q| ~ Wqo- Graphene sheet incorporates 
80000 atoms at 50 K. 



sheets. 

To compute the free-energy, we employed a well-known 
method (c.f., Haile [15]) in which one defines a continu- 
ous variable A for distinguishing two different states [H] . 
Suppose that by varying an external parameter, such as 
slow compression and relaxation of the graphene, the sys- 
tem can go from an initial state i (rough) to a final state 
/ (rippled) . When the parameters are changed infinitely 
slowly along some path from i to /in the parameter space, 
then the total work W performed on the system is equal 
to the Helmholtz free-energy difference between the ini- 
tial and final configurations. In contrast, when the pa- 
rameters are switched along the path at a finite rate, 
Jarzynski found that [T^: 



AA 



^lnexp(-^IF) 



(1) 



where overbar denotes an average over an ensemble of 
measurements of W. We ran the MD simulation to very 
long times in order to slowly pass the intermediate qua- 
sistable states. In practice, for every step of compression- 
relaxation, we checked whether the system was in equi- 
librium. We ensured the existence of the true equilibrium 
condition by checking the stability of the internal-energy 
fluctuations. Eventually, the problem of calculating AA 
is the same as calculating the averaged W. 

In Fig. 2 the free-energy differences ^ripple — trough is 
given as a function of T. We used the numerical results 
for a graphene sheet incorporating N = 80000 atoms 
(400 X 200) at various temperatures. It appears that the 
ripples are stable at low temperatures, namely, below 
Tc ~ 90 K, such that above T^ the rough state is more 
stable. This feature is in good agreement with recent 
experimental observation [lOj . Here, we would like to 
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FIG. 2: (Color online) The dependence of the free-energy and 
internal-energy (entropy) differences as a function of temperature. 
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FIG. 3; (Color online) Probability distribution function of 
the total energy E for the rippled and rough states for T = 
30,50,70,85,100 and 300 K (from top to bottom). The dashed 
curves are the Gaussian probability distribution function. For clar- 
ity, the PDFs were shifted upward. 



point out that the potential energy of the carbon-carbon 
interaction in the compression process is different from 
that in the relaxation process since the relative positions 
of the atoms in these two configurations are different. 
Note that the morphology of the surface depends strongly 
to the potential energy. 

We also tested the dependence of Tc on the size of 
the samples by simulating the systems with 600 x 200 
and 800 x 200 atoms, and found their characteristic tem- 
peratures to be Tc ~ 115 and Tc ~ 140, respectively. 
Moreover, we found that the wavelength of the ripples 
depends on T as A ~ 35 ln(T) — 55, but it does not de- 
pend on the system size. Accordingly, we calculated the 



entropy difference of the two states and showed that for 
temperatures less than Tc, the rippled state has a higher 
entropy and is stable. We plotted the internal energy 
difference of the two states is shown in Fig. 2. The value 
of Tc can be also determined from the local stored stress 
on a graphene sheet; we will report the results elsewhere 

m- 

We may expect that similar to second-order phase 
transitions the probability distribution function (PDF) 
of the total internal energy possesses different shapes for 
rough and ripple states, and exhibit non-Gaussian be- 
havior. In Fig. 3 the PDF of the total internal en- 
ergy E for the ripple and rough states are presented 
for T = 30,50,70,85,100 and 300 K. To calculate the 
PDF, we used 200 ensembles of roughened and rippled 
graphenes, incorporating 400 x 200 atoms. We observed 
that the PDF has a Gaussian form for both states indicat- 
ing that there is no longer second-order phase transition 
in the system. We also checked the Gaussian nature of 
the PDF by using the x^ test [5T]. The dashed curves 
represent the Gaussian PDF. 

As we argued earlier, there are at least two states for 
graphene sheets for a given temperature. A question 
raised is, whether or not, there is any possibility of a 
transition from one state to another? One possible way 
for such a transition with fixed graphene sheet size is to 
increase the temperature. For this purpose we simulated 
the graphene sheet with 80 x 40 atoms and, after carrying 
out the compression-relaxation process, the ripple struc- 
ture appeared at T = 55 K (see the upper figure of Fig. 
4). We then increased the temperature very slowly. As 
shown in Fig. 4, the ripples begin to disappear at high 
temperatures. At T = 55 K, we have almost two wave- 
length of the ripples; however, at higher temperatures 
there is one wavelength at T = 320 K, and finally one 
half of the wavelength at T = 520 K remains. The final 
step may be called rough state. Accordingly, the energy 
barrier of two states may be estimated by 465fcB = 0.04 
eV (or — 0.0125 meV per particle) for a sample with 3200 
atoms, where kg is the Boltzmann's constant. Such a 
transition has been observed experimentally in [10] . They 
argued that the disappearing of ripples in high tempera- 
ture is due to the fact that graphene has negative thermal 
expansion coefhcient. 

As mentioned above, for temperatures less than Tc, the 
free-energy of the rippled state is smaller than that of the 
free-energy of the roughened graphene sheet. However, 
there is a possibility of having a transition from the rip- 
pled state to the roughened state, due to a tunneling-type 
transition. To detect this transition, we checked that the 
height fluctuations variance of the rippled graphene sheet 
is stable with time. A sample size of 80 x 40 atoms was 
used at T = 55 K. The simulations showed that there 
is no transition from the rippled to roughened state at a 
constant temperature T less than Tc, at least up to avail- 
able time scales in the MD simulations. The probability 



for such transition is about exp(— 465/r) for a sample 
with size 3200 atoms. We remind that the variance of 
height fluctuations in graphene are about 5 A and 1 nm, 
for rough and rippled states, respectively. 
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FIG. 4: (Color online) Transition from rippled state (upper snap- 
shot) to the rough state, due to the increasing of the temperature 
from 55K to 520K. 



In summary, we have used a compression-relaxation 
mechanism to produce rippled structures on graphene 
sheets. The constructed ripples survive even though the 
system is relaxed to its initial position. In the closed-path 
loop, we calculated the total work and, hence, the free- 
energy difference of the rippled and roughened states. 
Our numerical results show that for sample with 400 x 200 
atoms and below Tc ~ 90 K, the rippled surface is sta- 
ble and the entropy of the ripples should be larger than 
that of the rough state. However, above Tc the rough 
state is more stable. The rippled and rough structures 
are also related to the morphology of such systems and 
we, therefore, expect that the our simulations yield the 
correct and new results for the free-energy of the rippled 
and roughened graphene. We have done similar simu- 
lations for a bilayer graphene and observed that, for a 
given temperature, the wavelength of the static ripples 
are larger than that for a monolayer graphene. We will 
report the results for bilayer graphene elsewhere. 
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